Optimal detection of losses by thermal probes 
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We consider the discrimination of lossy bosonic channels and focus to the case when one of the 
values for the loss parameter is zero, i.e., we address the detection of a possible loss against the 
alternative hypothesis of an ideal lossless channel. This discrimination is performed by inputting 
one-mode or two-mode squeezed thermal states with fixed total energy. By optimizing over this 
class of states, we find that the optimal inputs are pure, thus corresponding to single- and two-mode 
squeezed vacuum states. In particular, we show that for any value of the damping rate smaller than 
a critical value there is a threshold on the energy that makes the two-mode squeezed vacuum state 
more convenient than the corresponding single-mode state, whereas for damping larger than this 
critical value two-mode squeezed vacua are always better. We then consider the discrimination in 
realistic conditions, where it is unlikely to have pure squeezing. Thus by fixing both input energy 
and squeezing, we show that two-mode squeezed thermal states are always better than their single- 
mode counterpart when all the thermal photons are directed into the dissipative channel. Besides, 
this result also holds approximately for unbalanced distribution of the thermal photons. Finally, we 
also investigate the role of correlations in the improvement of detection. For fixed input squeezing 
(single-mode or two- mode) , we find that the reduction of the quantum Chernoff bound is a monotone 
function of the two-mode entanglement as well as the quantum mutual information and the quantum 
discord. We thus verify that employing squeezing in the form of correlations (quantum or classical) 
is always a resource for loss detection whenever squeezed thermal states are taken as input. 

PACS numbers: 03.67.-a, 42.50.Dv, 42.50.Ex 



I. INTRODUCTION 

One of the main obstacles to the development of quan- 
tum technologies is the decoherence associated to losses 
and absorption processes occurring during the propaga- 
tion of a quantum signal. The description of the dynam- 
ics of systems subject to noisy environments as well 
the detection, quantification and estimation of losses and, 
more generally, the characterization of lossy channels at 
the quantum level, received much attention in the recent 
years [1-4]. An efficient characterization of decoherence 
is relevant for quantum repeaters Q , quantum memories 
QED systems [TJ, superconducting quantum 



quantum teleportation [9j , quantum cryptog- 



|6|J, cavity 
circuits |8|| . qu ; 

raphy [l(J EH and secret- key capacities [IH, [l3[ . 

In this paper we address the discrimination of lossy 
channels, i.e., we consider a situation where the loss 
(damping) rate of a channel may assume only two pos- 
sible values and we want to discriminate between them 
by probing the channel with a given class of states. In 
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particular, we address the discrimination of lossy chan- 
nels for bosonic systems using squeezed thermal states 
as probing states, and focus attention to the case when 
one of the values for the loss parameter is zero, i.e., we 
address the detection of a possible loss against the al- 
ternative hypothesis of an ideal lossless channel. Such 
kind of discrimination is crucial since a recent analysis 
14j has revealed the importance of assessing the devia- 
tion from ideal conditions, i.e., the identity channel, in 
implementing large-scale quantum communication. 

Despite the discrimination of lossy channels has been 
already considered in the literature, the approach of this 
paper is novel. Previous works have in fact considered 
this kind of discrimination by constraining the energy ir- 
radiated over the unknown lossy channel but not the total 
energy employed in preparing the input states. For in- 
stance, in the quantum reading of digital memories [l5| . 
the discrimination of lossy channels has been analyzed 
by fixing the mean total number of photons irradiated 
over the channel, independently on the number of prob- 
ing modes (this approach has been also followed by the 
recent Ref. fl(| ) . Previously, in the quantum illumina- 
tion of targets [17l [Tc| . the channel discrimination was 
performed by fixing the mean number of photons in each 
of the modes probing the channel. In both these models 
there was no restriction for the energy involved in the 
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use of ancillary modes. Our approach considers the dis- 
crimination of lossy channels by constraining the total 
energy of the input state, thus including both the prob- 
ing mode (irradiated over the unknown lossy channel) 
and a possible ancillary mode (bypassing the lossy chan- 
nel and detected by the output measurement). Thus, 
while previous models were more focussed on restricting 
the energy irradiated over the channel, we address the 
problem from the point of view of the input source, i.e., 
considering the global effort in preparing this source. By 
fixing the input energy, we then optimize over an impor- 
tant class of Gaussian states, i.e., single- and two-mode 
squeezed thermal states. The choice of these states re- 
lies in their experimental accessibility, being routinely 
generated in today's quantum optics laboratories where 
they can be reliably controlled [1_9|. Furthermore, be- 
cause of the squeezing, they represent important exam- 
ples of non-classical states, i.e., states with non-positive 
P representation [2(j. This is another feature which di- 
versifies our study from previous works [ill 03 > where 
the main goal was the comparison between non-classical 
states and classical states (i.e., with positive P represen- 
tation). Note also that our problem, involving a discrete 
channel discrimination, is completely different from prob- 
lems of parameter estimation, where one has to infer a 
parameter taking continuous values. The estimation of 
the damping constant of a bosonic channel has been re- 
cently addressed by Refs. [2ll - l23 |. 

It is clear that, for a given input state, any problem of 
channel discrimination collapses into a problem of state 
discrimination [25H29J, where we have to compute the 
minimum error probability in identifying one of two pos- 
sible output states. By assuming M identical copies of 
the input state and a memoryless quantum channel, we 
have M output states which are exact replicas of two 
equiprobable states. In this case, the minimum error 
probability is well-approximated by the quantum Cher- 
noff bound (QCB) |30l - [34j |. Now the crucial step is to 
vary the input state, trying to optimize the value of the 
output QCB. In the case of two lossy bosonic channels, 
this kind of optimization must be constrained, meaning 
that we have to fix some crucial parameters of the in- 
put states, in particular, their energy. As we have al- 
ready mentioned before, in our investigation we optimize 
the QCB on the class of single- and two-mode squeezed 
thermal states by fixing their total energy. Single-mode 
thermal states are sent through the lossy channel, while 
two-mode thermal states probe the channel with the one 
of the modes (probing mode) while bypassing the chan- 
nel and assisting the output measurement with the other 
mode (reference mode). In this scenario, we find that 
the pure version of these states, i.e., single- and two- 
mode squeezed vacuum states, are optimal for detecting 
losses, i.e., discriminating a lossy from a lossless channel. 
Furthermore, we are able to show that, for any value of 
the damping rate smaller than a certain critical value, 
there is a threshold on the energy that makes the two- 
mode squeezed vacuum state more convenient than the 



corresponding single-mode state. More interestingly, for 
damping rates larger than the critical value, the two- 
mode squeezed vacuum state performs always better than 
the single-mode squeezed vacuum state with exactly the 
same energy. 

In order to stay close to schemes which are feasible 
with current technology, we also analyze the effect of the 
mixedness in the probe states. In this case, we study the 
channel discrimination by fixing not only the total energy 
of the input state but also its total amount of squeez- 
ing. Then, we are able to show that two-mode squeezed 
thermal states are always better than their single-mode 
counterpart when all the thermal photons are directed 
into the dissipative channel. We have numerically verified 
that this result also holds, approximately, for unbalanced 
distribution of the thermal photons. Finally, we also in- 
vestigate the role of correlations in the improvement of 
loss detection. In order to quantify correlations, besides 
entanglement and mutual information, we exploit the re- 
cent results on the quantum discord |35l44G |. which has 
been defined with the aim of capturing quantum correla- 
tions in mixed separable states that are not quantified by 
entanglement. Thus, at fixed input squeezing, we study 
the reduction of the QCB as a function of various corre- 
lation quantifiers, i.e., quantum mutual information, en- 
tanglement and quantum discord. This analysis allows us 
to conclude that employing squeezing in the form of cor- 
relations, either quantum or classical, is beneficial for the 
task of loss detection whenever squeezed thermal states 
are considered as input probes. 

The paper is structured as follows. In Sec. |TT] we re- 
view the discrimination of quantum states together with 
the general definition of the QCB. In Sec. IHII we give the 
basic notation for Gaussian states and in Sec. |TV] we re- 
view the formula of the QCB for Gaussian states. Then, 
in Sec. |V] we discuss the discrimination of lossy channels 
by considering single-mode and two-mode squeezed ther- 
mal states. In particular, we provide the details on how 
to compute the QCB for distinguishing an ideal lossless 
channel from a lossy channel. Section I VII reports the 
main results regarding the single- and two-mode states 
as a function of the total energy and squeezing. Finally, 
in Section I VIII we analyze the role of correlations in en- 
hancing the discrimination. Section rVIIII closes the paper 
with some concluding remarks. 



II. QUANTUM STATE DISCRIMINATION AND 
QUANTUM CHERNOFF BOUND 

As we have mentioned before, the problem of quantum 
channel discrimination collapses to the problem of quan- 
tum state discrimination when we fix the input state. As 
a result, mathematical tools such as the quantum fidelity 
and the QCB are fundamental in our analysis. Despite 
it has been introduced only very recently, the QCB has 
been already a crucial tool in several areas of quantum 
information: it has been exploited as a distinguishabil- 
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ity measure between qubits and single-mode Gaussian 
states [lil Hl[ , to evaluate the degree of nonclassicality 
for single-mode Gaussian states [4l[ or the polarization 
of a two-mode state [42| ■ It has also been applied in the 
theory of quantum phase transitions to distinguish be- 
tween different phases of the XY model at finite temper- 
ature [43d , and to the discrimination of two ground states 
or two thermal states in the quantum Ising model [44| . 
For continuous variables systems the quantum discrim- 
ination of Gaussian states is a central point in view of 
their experimental accessibility and their relatively sim- 
ple mathematical description (4|||46j]. In fact, in the case 
of Gaussian states, the QCB can be computed from their 
first and the second statistical moments. A first formula, 
valid for single-mode Gaussian states, was derived in [30l |. 
Later, Ref. [34[ provided a general closed formula for mul- 
timode Gaussian states, relating the QCB bound to their 
symplectic spectra. Furthermore, from these spectra one 
can derive larger upper bounds which are easier to com- 
pute than the QCB [Hj]. 

In this section, we start by establishing notation and 
reviewing the problem of quantum state discrimination, 
together with the general definition of QCB. Then, from 
the next section, we will specialize our attention to the 
case of Gaussian states and we will review the formula 
for Gaussian states in Sec. IIV1 

In its simplest formulation, the problem of quantum 
state discrimination consists in distinguishing between 
two possible states, pa and ps, which are equiprobable 
for a quantum system. We suppose that M identical 
copies of the quantum system are available. Then, we 
have the following equiprobable hypotheses, Ha and Hb, 
about the global state 



H A : pf = p A 



PA 



H B : Pb 1 = P b 



PB ■ 



M 

In order to discriminate between these two hypothe- 
ses, one can measure the global system by using a 
two-outcome positive operator valued measure (POVM) 
{E A , E B }, with Ea + Eb =1 and E A ,E B > 0. After ob- 
serving the outcome j = A or B, the observer infers that 
the state of the system was pj 1 . The error probability 
of inferring the state p^ when the actual state is p^ is 
thus given by the Born rule Pjk = Tr [pj^.Ej] . As a re- 
sult, the optimal POVM for this discrimination problem 
is the one minimizing the overall probability of misiden- 
tification, i.e., P e = \{Pba + Pab)- Since Ea = I — E B , 
we have 

P e = ±Tr\p%E B ] + \t±\p%E a ] 



where 



= ~ (1 - Tr [E B A]) 



A = Pb I -Pa 1 



(1) 



is known as the Helstrom matrix |25| . Now, the er- 
ror probability P e has to be minimized over Eb- Since 
Tr[A] = 0, the Helstrom matrix has both positive and 
negative eigenvalues and the minimum P e is attained if 
Eb is chosen as the projector over A + , i.e., the positive 
subspace of A. Assuming this optimal operator we have 
Tt[E b A] = Tr[A+] = ±Tr|A| with |A| = VAU. Thus 
the minimal error probability is given by 

1 



where 



P* = $[l-T(j%,p%j\ 



T(p,<T) = -Tr\p-a\ 



is the so-called trace distance. The computation of the 
trace distance may be rather difficult. For this reason, 
one can resort to the QCB that gives an upper bound to 
the probability of error P e [30l - l34| 



P, < 



Q 



M 



where 



(2) 



(3) 



The bound of Eq. ([2]) is attainable asymptotically in the 
limit M 00 as follows from the results in [3l], H2 . One 
may think that the trace distance has a more natural 
operational meaning than the QCB. In spite of this, it 
does not adapt to the case of many copies; indeed, one 
can find states p, a, p' , a' such that 

T{p, a) < T(p', a') but T(p' M , a' M ) < T(p M , a M ) . 

By contrast, the QCB does resolve this problem since 

Q{p,a)<Q{p',a') => Q(p M ,a M ) < Q(p' M ,<r' M ) . 

Because of this property, the minimization of the QCB 
over single-copy states (p and a) implies the minimization 
over multi-copy states (p M and a M ). This is true as long 
as the minimization is unconstrained or if the constraints 
regard single-copy observables (e.g., the mean energy per 
copy). 

Finally, note that there is a close relation between the 
QCB and the Uhlmann fidelity 



F(pa,Pb) = Tr 



PAPB^PA 



which is one of the most popular measures of distin- 
guishability for quantum states. In fact, for the single- 
copy state discrimination (M = 1) we have [30ll33 . l47ll48j 



l-y/l-F(p A ,p B ) 



Q 

< P e < 7T < 
- - 2 ~ 



\/F{pa,Pb) 



(4) 
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More generally, by exploiting the multiplicativity of the 
fidelity under tensor products of density operators, i.e., 

F(p A ® cr Al p B ® ctb) = F(p A ,p B )F(a A ,a B ) , 

we can write 

F{pf, PB *) = F{p A ,p B ) M = F M . 

This leads to the general multi-copy version of Eq. (HJ 
which is given by [49( 



1 - VI - F M Q M F I 
1 < p < ^_ < 

2 - 2 ~ 2 



(5) 



From the previous inequalities, it is clear that the QCB 
gives a tighter bound than the quantum fidelity. How- 
ever, if one of the two states is pure, then the QCB just 
equals the fidelity, i.e., we have 

Q(pa,Pb) = F{p A ,p B ) = Tr[p A p B ] . 



III. GAUSSIAN STATES 

In this section we give the basic notions on bosonic 
systems and Gaussian states, ending with the definition 
of squeezed thermal states. 

An n-mode bosonic system is described by a tensor- 
product Hilbert space R® n and a vector of canonical 
operators R = (qi,pt, ■ ■ ■ , q n ,Pn) T satisfying the com- 
mutation relations 



where I, m = 1 , • 
symplectic form 



[Ri,R m ] = iQi m , 

, 2n and fi/ m are the elements of the 



k=l 



1 

-1 



(6) 



Alternatively, we can use the mode operators which 
are given by the cartesian decomposition of the canonical 
operators, i.e., 



1 



1 



These operators satisfy the commutation relations 
[a k ,a\,] = 5 k k' with k, k' = 1, • • • , n. 

An arbitrary quantum state p of the system is equiva- 
lently described by the characteristic function 

x[p](X)=Tr[ P D(X)} 

where D(X) — <8>£ =1 -Dfc(Afc) is the n-mode displace- 
ment operator, with A = (Ai, . . . , A„) T , A^ £ C, and 
-Dfc(Afe) = expjAfeaJ. — A^afe} is the single mode displace- 
ment operator. A state p is called Gaussian if the corre- 
sponding characteristic function is Gaussian 



\[p](A\ = cxi)<| -iA T crA- 



X 1 CIA 



(7) 



where A is the real vector 

A = (ReA l5 ImAi, . . . , ReA„, ImA„) T . 

In this case, the state is described by its first two statis- 
tical moments, i.e., the vector of mean values X and the 
covariance matrix (CM) cr, whose elements are defined 
as 



X t =(R) 

Vim =^\{Rl)Rm\) — (Rl)(Rr, 



(8) 



where {A, B} = AB + BA denotes the anti-commutator, 
and (O) = Tr[pO] is the mean value of the operator O. 

In the remainder of this section, we consider only zero- 
mean Gaussian states, i.e., Gaussian states with X = 0, 
which are therefore fully specified by their CM. The prop- 
erties of these states may be expressed in very simple 
terms by introducing the symplectic transformations. A 
matrix S is called symplectic when preserves the sym- 
plectic form of Eq. ©, i.e., 



sns J 



n. 



Then, according to the Williamson's theorem, for every 
CM <t, there exists a symplectic matrix S such that 



a = SWS T 



(9) 



where 



k=l 



1 
1 



and the dk's are called the symplectic eigenvalues of cr. 
The physical statement implied by the decomposition of 
Eq. ([9]) is that every zero-mean Gaussian state p can be 
obtained from a thermal state by performing the unitary 
transformation Us associated with the symplectic matrix 
S, i.e., we have 



where v — v\® 
thermal states 

Vk = 



p = U s uUl 
) v n is a tensor product of single- mode 



n k 



nk + 1 



\m)k(m\ 



with average number of photons given by fik = dk — 1/2. 
For a single-mode system the most general zero-mean 
Gaussian state may be written as 

where S(() — exp{i(£a^ 2 — C* fl2 )} is the single-mode 
squeezing operator and C = r ^ G C. The corresponding 
covariance matrix is given by 



a c 
c b 



(10) 
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where 



a = (n H — ) [cosh(2r) — sinh(2r) cos < 
b = (h + —) [cosh(2r) + sinh(2r) cos < 
c = (h + — ) sinli(2r) sin0 . 



(11) 



In particular, we can consider the case of a real squeezing 
parameter, e.g., by fixing ( = —r [5(J. In this case, the 
previous expressions of Eq. (jlip simplify into the follow- 
ing 

a = ^(2«+ l)exp(2r) 

b = i(2n + l)exp(-2r) 

c = 0. (12) 

This state defines the single-mode squeezed thermal 
state. It depends on two real parameters only, i.e., we 
have 

p = S(r)vS' (r) = p(r, n) . 

In particular, for n — the state is pure and corre- 
sponds to a single-mode squeezed vacuum state p(r, 0) = 
S(r)|0)<0|St( r ). 

Now let us consider two-mode (zero-mean) Gaussian 
states. They are completely characterized by their 4x4 
CM 



A C 
C T B 



(13) 



where A, B and C are 2x2 blocks. It is useful to intro- 
duce the symplectic invariants 



h = det A , h = det B . 
^3 = det C , I4 — det er 



(14) 



By means of these invariants, we can simply write the 
two symplectic eigenvalues as 



d± = 



'A± 



where A = h + I2 + 2/3 [5l|, [52|. By means of local 
symplectic operations, the CM of Eq. (Tl~3")) can be recast 
in the standard form, where the three blocks A and B 
are proportional to the identity and C is diagonal [5lj ]. 
In the particular case of a two-mode squeezed thermal 
state, we can write 



1 ( Ah Ca z 



2 V Ca z BI 2 



(15) 



where 



A = cosh(2r) + 2ni cosh 2 r + 2n 2 sinh 2 r 

B = cosh(2r) + 2ni sinh 2 r + 2t%2 cosh 2 r 

C= (l + ni+n 2 )sinh2r, (16) 



with I2 the 2x2 identity matrix and cr z — diag(l, — 1) 
the z-Pauli matrix. This corresponds to considering a 
density operator of the form 

P = S 2 {r){v 1 ®V2)S 2 {r)\ 

where S>z(r) = exp{r(a 1f b^ — ab)} is the two-mode squeez- 
ing operator. This state depends on three real parame- 
ters: the squeezing parameter and the two thermal num- 
bers, i.e., we have 

p = p(r,ni,n 2 ) ■ 

In particular, for n\ = n% — the state is pure and corre- 
sponds to a two-mode squeezed vacuum state p(r, 0, 0) = 
S2(r)(\0)i <0| ® |0> a <0|)4(r). 



IV. QUANTUM CHERNOFF BOUND FOR 
GAUSSIAN STATES 



Here we review the formula of the QCB for multimode 
Gaussian states 34]. In particular, we adapt this formula 
to our notation and physical units (here the vacuum noise 
is 1/2, while in Ref. [34j it was equal to 1). Let us consider 
two Gaussian states p (with statistical moments X and 
<t) and p' (with statistical moments X and er'). The 
CMs of these two states can be decomposed as 

■ 71 n 



<T = SW(n 1 ,---,n n )S 1 (17) 
*' = S' W(h[,--- ,n' n ) S' T , (18) 
where {fik} and {n f k } are their thermal numbers, and 

n 

W{x u --- , Xn ) = 0(2x fe + l)I 2 . 

k=l 



Now let us define the functions 



and 



G s (x) 



A s (x) 



1 



(x + iy 



{x + iy 

Then, the QCB is given by 



Q 



inf Q s , 

0<s<l 



where 



Q s = 



n. 



: exp 



1 



VdetEl V 2 
In the formula of Eq. P^|) . we have d = X — X', 



(19) 



n s = J] G s {n k )G^ s (n' k 



k) j 



fc=l 



and 



S s =S W[A s (ni), ■ ■ ■ ,A s (n n )} S T 

+ S' W[A^ s (n[),--- ,Ai- s «)] S' T . 
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A. Discrimination of squeezed thermal states 

For the discrimination of squeezed thermal states, the 
previous formula simplifies a lot. First of all, since they 
are zero-mean Gaussian states, we have d = and, there- 
fore, the exponential factor in Eq. (|T9l) disappears. Then, 
the symplectic decompositions in Eqs. (fTT]) and (|18[) arc 
achieved using symplectic matrices S and S' which are 
just one-parameter squeezing matrices, i.e., S = S(r) 
and S' = S'(r'). 

Thus, let us consider the discrimination of single- mode 
squeezed thermal states p = p(r,n) and p' = p'(r',n'). 
In this case, the QCB can be computed using 



where 



and 



Qg= n*(n,n') 

y/dcf£ s (r,n,r',n') ' 

U s (n,n') = G s (n)Gi- s {n') , 

•£ s (n,n',r,r') = S(r)W[A s (n)]S(r) 

+ S(r')W r [A 1 _ s (n')]S(r 



(20) 



For the discrimination of two-mode squeezed thermal 
states p — p(r, n\,n2) and p' — p'(r', n[, n' 2 ), we can use 



Qa 



U s (ni,n2,n' v n' 2 ) 



(21) 



where 

n s (ni,n 2 ,n' 1 ,n' 2 ) = G s (ni)G s (n2)Gi- s (n' 1 )Gi- s (n 2 ) , 
and 

E^r.ni.nay.ni.ni) = S(r)W[A s (n 1 ),A s (n 2 )]S(r) T 
+S(r')W[A 1 _ s (n' 1 ),A 1 _ s (n 2 )]S(r') T . 



V. DETECTION OF LOSSES BY THERMAL 
PROBES 

In what follows, we study the evolution of a Gaus- 
sian state in a dissipative channel £r characterized by 
a damping rate T, which may result from the interac- 
tion of the system with an external environment, as for 
example a bath of oscillators, or from an absorption pro- 
cess. We consider the problem of detecting whether or 
not the dissipation dynamics occurred. Given an input 
state p, this corresponds to discriminating between an 
output state identical to the input p, and another out- 
put state storing the presence of loss £r(p)- Lossy chan- 
nels are Gaussian channels, meaning that they tansform 
Gaussian states into Gaussian states. Furthermore, if 
the input is a squeezed thermal state, then the output 




FIG. 1: Single- and two-mode schemes for the detection of 
losses. Top: a single-mode squeezed thermal state p enters the 
lossy channel with damping rate F. A measurement apparatus 
detects the output state p' . Bottom: the lossy channel acts on 
the probing mode (1) of a two-mode squeezed thermal state 
p, while the reference mode (2) bypasses the channel. The 
output state p' of both the modes is then measured. 



state is still squeezed thermal (this is discussed in detail 
afterwards). 

In general, we consider the schematic diagram depicted 
in Fig[TJ In order to detect loss, we consider either a 
single-mode squeezed thermal state evolving in the lossy 
channel with parameter T followed by a measurement at 
the output, or a two-mode squeezed thermal state with 
the damping process occurring in only one of the two 
modes (the probing mode), followed by a measurement 
on both of the modes. Our aim is to minimize the error 
probability in discriminating between the ideal case T = 
and the lossy case T > 0. In the next section, this 
will be done by fixing some important parameters of the 
input state, such as total energy and squeezing [53|. 

The propagation of a mode of radiation in a lossy chan- 
nel corresponds to the coupling of the mode o with a zero 
temperature reservoir made of large number of external 
modes. By assuming a Markovian reservoir and weak 
coupling between the system and the reservoir the dy- 
namics of the system is described by the Lindblad Master 
equation [H| 



p = -C[a]p 



(22) 



where £[a]p = lapa) — a^ap — pa^a. The general solution 
may be expressed by using the operator-sum representa- 
tion of the associated completely-positive map, i.e., upon 
writing r\ — e~ rt , we have 



vi 



where 



{l-r]) m m i (n t n _ m) 
■ / 



and q is the initial state. 
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A. Single-mode case 

Let us now start with single-mode states. Eq. (f22j) can 
be recast into a Fokker-Planck equation for the Wigner 
function in terms of the quadrature variables q and p. 



w = - + 0£<roo0£] w 



where X = (q,p) T , dx = (d q ,d p ) T and we introduced 
the diffusion matrix a x = diag(l/2, 1/2). Solving the 
equation for the Wigner function of a single-mode Gaus- 
sian state one can obtain the evolution equation for the 
CM cr. This is given by [H 



-r( CT 



>) 



which yields to 
cr(t) 



-rt 



(T + (1 - e 



-rt\ 



The latter equation describes the evolution of an initial 
Gaussian state with CM cr towards the stationary state 
given by the Gaussian state of the environment with CM 
(Too. For simplicity, from now on we omit the index of 
(To, we replace ait) — > cr' , and we insert the time t into 
the damping parameter T. Thus, the evolved CM of the 
single mode case simply reads 

cr' = e~ T a + (1 — e~ r ) cr oo . 

Now let us consider the specific case of an input 
squeezed thermal state p — p(r, tit) with squeezing r and 
thermal number nx- According to Eqs. (fTU)) and (|12p . its 
CM is given by 



a 
b 



with 

a=^(l + 2n T )e 2r , b = ± (1 + 2n T ) e' 2r . 
At the output of the channel, the state p' has CM 



(23) 



Now, the discrimination between a lossless (T = 0) and a 
lossy channel (r > 0) corresponds to the discrimination 
between the input state p = p(r, ut) and the output one 
p' = p'(rr,nr). In order to estimate the error proba- 
bility affecting this discrimination, we can compute the 
quantum Chernoff bound. This is achieved by replacing 

(r, n) — > (r, riy) and (r', n') — > (tt, nr) 

in Eq. (|20]>. 

B. Two-mode case 

According to the scheme of Fig. Q] the map describing 
the evolution of a two-mode state is £ r <£> I, where the 
lossy channel £r acts on the probing mode while the iden- 
tity channel 1 acts on the reference mode. At the level 
of the CM it corresponds to the following transformation 

cr' = (e- r / 2 I 2 ffi I 2 ) cr (e- r/2 I 2 © I 2 ) 

+ {I 4 -e- r h®h)cr 00 . (27) 

As input state, let us consider a two-mode squeezed ther- 
mal state p = p(r, nx 1 , n^). Its CM is provided in 
Eq. (|15p with the elements given in Eq. (|16[) by replacing 
(r,ni,n 2 ) — > (r, , i-t 2 ) ■ The CM of the output state 
can be derived using the Eq. (|27|) . This CM can be put 
in the normal form of Eq. (1151) with elements given by 
Eq. p^|) by replacing (r,ni,n 2 ) — > (rp, uy 1 , ^r 2 )- Here 
the squeezing parameter rr and the thermal numbers nr 1 
and nr 2 are function of the input parameters r, nr 1 and 
n,T 2 (the explicit expression is too long to be shown here) . 
Thus, the output state is still a two-mode squeezed ther- 
mal state p' — p'(rr,ni' 1 ,nr 2 ). 

As before, the discrimination between a lossless (r = 
0) and a lossy channel (r > 0) corresponds to the discrim- 
ination between the input state p = p(r, , tit 2 ) and 
the output one p' = p'(rr, nr 1 , nr 2 ). The error probabil- 
ity affecting this discrimination is estimated by the QCB 
which is computed by replacing (r, n%, n 2 ) — > (r, , tit 2 ) 
and (r'jn^n^) — > (rr,nr 1 ,nr 2 ) in Eq. (|2"TT) . 



a r 
b r 



where 



a r = \ (1 + 2n r ) e 2rr , b r = ± (1 + 2n r ) e" 2 ^ 



and 



n r = ^dct[cr'} - 1/2 



r r = - log 



e- r a + (l-e- r )/2 
e- r fe+(l-e- r )/2 



(24) 

(25) 
(26) 



Thus, we still have a squeezed thermal state p' — 
p'(rr,nr) with squeezing rr and thermal number nr- 



VI. OPTIMIZATION OF THE THERMAL 
PROBES 

In this section, we optimize the discrimination of a loss- 
less (r = 0) from a lossy channel (r > 0) by maximizing 
over thermal probes, i.e., single- and two- mode squeezed 
thermal states. For this sake, we evaluate the QCB as 
a function of the most important parameters of the in- 
put state, i.e., its total energy and squeezing. In our 
first analysis, we show that for fixed total energy, single- 
and two-mode squeezed vacuum states are optimal. In 
particular, we show the conditions where the two-mode 
state outperforms the single-mode counterpart. Then, by 
fixing both the total energy and squeezing, we will find 
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the optimal squeezed thermal state. According to Sec. |TT] 
the minimization of the QCB over single-copy states im- 
plies the minimization over multi-copy states (when the 
minimization is unconstrained or subject to single-copy 
constraints). This implies that finding the optimal in- 
put state p at fixed energy automatically assures that 
p ® p ® ■ ■ ■ is the optimal multi-copy state at fixed en- 
ergy per copy when we consider a multiple access to the 
unknown (memoryless) channel. 

In order to perform our investigation we introduce a 
suitable parametrization of the input energy. Given a 
single-mode squeezed thermal state p = p(r, tit), its en- 
ergy (mean total number of photons) can be written as 



Ni = tit + ns + 2nsnT 



(28) 



where nx accounts for the mean number of thermal pho- 
tons, ns = sinh 2 r quantifies the squeezing, and nsnr is 
a cross term. Alternatively, we can introduce a squeezing 
factor j3i € [0, 1] such that 



n s = p x Nx 

(1 - PM 
1 + 2/3iJVi 



nT 



(29) 
(30) 



Thus the single-mode squeezed thermal state can be 
parametrized as p — p(Ni,j3{), i.e., in terms of its to- 
tal energy Ni and the squeezing factor /3i- Note that 
for fix — the state is completely thermal with energy 
Ni = n T , while for /3x = 1 the state is a squeezed vacuum 
with energy N\ = ns- In our problem of loss detection 
(r = versus T > 0), we denote by Qi(iVi, the out- 
put QCB which is computed by using the input state 
p = p(N u pi). 

Now, given a two-mode squeezed thermal state p — 
p(r, riTi , tit 2 ) , its total energy can be written as 



N 2 = n Tl + n T2 + 2n s + 2n s (n Tl + n T2 ) , 



(31) 



where {nr 2 ) quantifies the thermal photons in the 
probing (reference) mode, ns = sinh 2 r quantifies the 
two-mode squeezing energy, and the last energetic term is 
a cross term. In this case, besides the squeezing factor (3 2 , 
we can also introduce an asymmetry parameter 7 G [0, 1] 
which quantifies the fraction of thermal energy used for 
the probing mode. In other words, we can write 



f 



ns 



n Tl = 7 - 



n T2 



-J2N2 

( 1 - h)N 2 
1 + P2N2 

(1 - h)N 2 



(1-7) 



1 + P2N2 



(32) 
(33) 
(34) 



Thus the two-mode squeezed thermal state can be 
parametrized as p — p(N 2 , fc, 7), i.e., in terms of the 
total energy N 2 , the squeezing factor /?2, and the asym- 
metry parameter 7. Note that for /3 2 = we have two 
thermal states, one describing the probing mode with 



thermal energy nji = 7 N 2 , and the other one describing 
the reference mode with thermal energy n? 2 = (l — "f)N 2 . 
For j3 2 = 1 we have instead a two-mode squeezed vacuum 
state with total energy N 2 = 2ns- In this case, the ther- 
mal energy is zero and 7 can be therefore arbitrary. In 
our problem of loss detection (r — versus T > 0), we 
denote by Q 2 (A r 2, 02, 7) the output QCB which is com- 
puted by using the input state p = p(A^2, 02, 7)- 



A. Optimal input at fixed total energy 

In our first investigation we fix the mean total number 
of photons of the input state. In other words we fix 



^1=^2=^. 



(35) 



Then we minimize the output QCB among single-mode 
and two-mode squeezed thermal states. As a first step 
we compute the optimal quantities 



Qi(N) :=mlQ x {N,Pi) 

Pi 

Q 2 (N) := inf Q 2 {N,p 2 ,l) 

02,1 

Then, we compare Qi(N) with Q 2 {N). 



(36) 
(37) 
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FIG. 2: Output QCB Qi(N,f5i) optimized over input single- 
mode squeezed thermal states p — p(N, /3i). From left to right 
we consider different values of the transmissivity: r\ = 0.1 (left 
panel), r\ = 0.5 (middle panel) and r\ — 0.9 (right panel). In 
each panel, we plot Qi(N,/3i) as function of the energy TV 
for different values of /9i. From top to bottom: /3i = 0.1 
(dashed line), /3i = 0.5 (dotted line) and Pi = 1 (solid line). 
The minimum curve is always achieved for /3i = 1, i.e., for an 
input single-mode squeezed vacuum state. 

According to our findings, in the Eqs. (131)1) and (|3T|) the 
infima are achieved for /3i = $2 = 1. This is numerically 
shown in Fig. [2] for the single- mode case and in Fig. [3] for 
the two-mode case. Thus, we have found that, at fixed 
input energy N , the optimal thermal probes are given by 
single- and two-mode squeezed vacuum states. In this 
case, the input state is pure and the QCB corresponds to 
the fidelity (which is the case when the s-overlap in Eq. 
(|3"|) is minimized for s approching the border). Let us 
adopt the transmissivity r\ = e" 
ing of the channel, so that T = 
sponds to 77 = 1, and T > (lossy channel) corresponds 
to < 77 < 1 . Then, for single-mode we can write 



r to quantify the damp- 
= (ideal channel) corre- 



Qi(N) = 



1 



^l + N(l-n 2 ) 



(38) 
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FIG. 3: Output QCB Q2 (N, fa , 7) optimized over input two- 
mode squeezed thermal states p = p(N,fa,j). From left to 
right we consider different values of the transmissivity: r\ — 
0.1, 0.5 and 0.9. From top to bottom, we consider different 
values of the asymmetry parameter 7 = 0, 0.5 and 1. In 
each panel, we then plot Q2 as function of the energy N for 
different values of fa. From top to bottom: fa = 0.1 (dashed 
line), fa = 0.5 (dotted line) and fa = 1 (solid line). The 
minimum curve is always achieved for fa = 1 corresponding 
to an input two-mode squeezed vacuum state. 



and for two-modes we derive 



Qi{N) = 



[2 + iV(l-0y)]' 



(39) 



In Fig. HI we show the behaviors of the single-mode QCB 
Qi(N) and two-mode QCB Q2(N) as function of the in- 
put energy N for several values of transmissivity r\ (or, 
equivalently, the damping rate T). As expected the dis- 
crimination improves by increasing the input energy N 
and decreasing the transmissivity 77. 

As we can see from Fig. [4j for a given value of the 
transmissivity 77, the two-mode QCB Q2(N) outperforms 
the single-mode QCB Qi(N) only after a threshold en- 
ergy. In fact, for any value of the transmissivity r\ 
larger than a critical value r\ c there is a threshold energy 
Nth — N t h {if) that makes the two- mode squeezed vacuum 
state more convenient than the single-mode counterpart. 
This threshold energy decreases for decreasing values of 
7/. In particular, for transmissivities less than the critical 
value r] c , the threshold energy becomes zero, i.e., the two- 
mode state is always better than single-mode state. We 
have numerically evaluated the critical value r] c ~ 0.296 
(corresponding to r c ~ 1.22). This phenomenon is fully 
illustrated in Fig. [5j where we have plotted the threshold 
energy as function of the transmissivity N t h = N t h{ii). 
For N > N t h (dark area), the optimal state is the two- 




10 15 20 



FIG. 4: (Color online). Single-mode QCB Qi(N) (solid lines) 
and two- mode QCB Qz(N) (dashed lines) as a function of 
the input energy N for different damping rates. From top 
to bottom r = 0.1,0.3,1 (red, green and blue, respectively) 
corresponding to r\ ~ 0.9,0.74,0.37. By comparing curves 
with the same color (fixed damping F), we can see that Q2(N) 
outperforms Qi(JV) only after a certain value of the input 
energy TV. 




FIG. 5: Threshold energy as a function of the transmissivity 
Nth ~ Nth(ri) (solid curve dividing the dark and the white ar- 
eas). The dark area indicates the values of the energy N for 
which the two-mode squeezed vacuum state is optimal. The 
white region indicates where the single-mode squeezed vac- 
uum state is optimal. The dashed line denotes the behavior 
of the threshold energy Nth close to the critical transmissivity 
r/ c ~ 0.296. 



mode squeezed vacuum state, while for N < Nth (white 
area) it is the single-mode squeezed vacuum state. In 
particular, note that Nth = at 77 = r\ c . Close to the 
critical transmissivity we have [56} 



Nth ^ 4(77 - 77 c )+ 5.5(?7 -Tic) 2 



(40) 



B. Optimal input at fixed energy and squeezing 

It should be said that, in realistic conditions, it is un- 
likely to have pure squeezing. For this reason, it is im- 
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portant to investigate the performances of the squeezed 
thermal states by fixing this physical parameter together 
with the total energy. Thus, in this section, we fix both 
the input energy and squeezing, i.e., we set 



Nt = N 2 = N, 
Pi=p2=0 (0 < (3 < 1) . 



(41) 



Then, we compare the single-mode squeezed thermal 
state p = p(N,j3) with the two-mode squeezed thermal 
states p = p(N,/3,~f) for various values of 7. In other 
words, we compare Qi(N,(3) and Q 2 (N, /3, 7). 

For fixed N and (3, we find that the minimum of 
Q 2 {N, (3, 7) is achieved for 7 = 1 (easy to check numeri- 
cally) . This means that two- mode discrimination is easier 
when all the thermal photons are sent through the lossy 
channel. In this case we find numerically that 

Q 2 (N,I3,1) <Qi(N,/3), 

for every values of the input parameters N and /3, and 
every value of damping rate T in the channel. In other 
words, at fixed energy and squeezing, there is a two-mode 
squeezed thermal state (the asymmetric one with 7 = 
1) able to outperform the single- mode squeezed thermal 
state in the detection of any loss. In order to quantify 
the improvement we introduce the QCB reduction 

AQ = Q 1 (N,0)-Q 2 (N,P,1). 

The more positive this quantity is, the more convenient is 
the use of the two-mode state instead of the single-mode 
one. In Fig. [5] we show the behavior of AQ as function 
of the input energy and squeezing for two different values 
of the damping. As one can see from the plot, the QCB 
reduction is always positive. Its value increases with the 
energy while reaching a maximum for intermediate values 
of the squeezing. By comparing the two panels of Fig. [6j 
we can also note that the QCB reduction increases for 
increasing damping T (i.e., decreasing transmissivity). 




FIG. 6: (Color online) Density plot of the QCB reduction AQ 
as function of the input energy N and the squeezing /3. The 
left plot is for T = 0.1 and the right one for T = 0.9. 

Thus, we have just shown that, for fixed values of N 
and j3, the asymmetric two-mode squeezed thermal state 
(7 = 1) is the optimal thermal probe in the detection 



of any loss T. Here we also show that this is approxi- 
mately true for 7 < 1. In other words, we show that the 
inequality Q 2 {N, /3, 7) < Q\(N,(3) is robust against fluc- 
tuations of 7 below the optimal value 7=1. This prop- 
erty is clearly important for practical implementations. 
To study this situation, let us consider the 7-dependent 
QCB reduction 



AQ^Q^N, P)-Q 2 {N,p, 7 ) 



(42) 



In Fig. [7] we have specified this quantity for different val- 
ues of the asymmetry parameter 7 (each panel refers to 
a different value of 7). Then, for every chosen 7, we have 
computed AQ 7 over a sample of 10 3 random values of 
N, (3, and T (in each panel). As one can see from the 
figure, the quantity AQ 7 is approximately positive also 
when 7 is quite different from the unity. 
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FIG. 7: (Color online) QCB reduction AQ 7 for different val- 
ues of 7 (top left 7 = 0.99, top right 7 = 0.9, bottom left 
7 = 0.8, bottom right 7 = 0.7). In each panel, AQ 7 is com- 
puted over a sample of 10 3 random values of N, (3, and T. 



VII. ANALYSIS OF THE CORRELATIONS 

Since two-mode squeezed thermal states are able to 
outperform the single-mode counterpart under several 
physical conditions, it is natural to investigate this im- 
provement directly in terms of the correlations of the in- 
put state. The quantification of the correlations is real- 
ized by using the entanglement, the quantum discord and 
the quantum mutual information. In order to quantify 
the degree of entanglement of a two-mode Gaussian state, 
we can use the logarithmic negativity. Let us consider a 
bipartite Gaussian state with CM given in Eq. (fT5|) . It is 
easy to derive the symplectic eigenvalues of the partially 
transposed state. These are given by 



d± = 



A - \/A 2 — 4/4 



where A = 1\ + I 2 — 2/ 3 , and the symplectic invariants 
Ii, I2, I3 and I4 are defined in Eq. (|T4")) . From the small- 
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est of these symplectic eigenvalues, we can compute the 
logarithmic negativity, which is equal to 

E = max{0,-log2d_} . 

A bipartite Gaussian state is entangled iff d_ < 1/2, so 
that the logarithmic negativity gives positive values for 
all the entangled states and otherwise. 

The quantum discord is defined as the mismatch of 
two different quantum analogues of classically equivalent 
expressions of the mutual information and may be used to 
quantify quantum correlations in mixed separable states. 
For a two-mode squeezed thermal state with CM as in 
Eq. (|15p . the quantum discord may be written as [38| 



D =h(y/I- 



h{d^) - h(d+) 



1 



(43) 



where 
h(x) = 



In.„ 



I+ 2 



x — — ) log \ x — — 



is the binary Shannon entropy. We have that, for < 
D < 1, the state may be either entangled or se par able, 
whereas all the states with D > 1 are entangled [3Sl . f39j j . 

Finally, the quantum mutual information, which quan- 
tifies the amount of total, classical plus quantum, corre- 
lations, is given by I = S(pa) + S(pb) — S(pab), where 
S(p) — — Tr[plogp] is the von Neumann entropy of the 
state p and pa(b) = ^b(A)[pab] are the partial traces 
over the two subsystems. For a two-mode squeezed ther- 
mal state with CM as in Eq. (|15[) the quantum mutual 
information can be computed using the formula 



I 



hwh) + hWh) - h{d+) - h{d-) 



For pure states the previous three measures are equiva- 
lent, whereas for mixed states, as in the case under inves- 
tigation in this section, they generally quantify different 
kind of correlations. Here we consider the QCB reduction 
AQj = Qi(N,f3) — Q%{N, /3, 7) between a single-mode 
squeezed thermal state p — p(N,j3) and a two-mode 
squeezed thermal state p — p(N,/3,^y) with 7 = 0.999. 
By fixing the input squeezing /3 and varying the input 
energy N, we study the behaviour of AQj as function 
of the three correlation quantifiers, i.e., quantum mutual 
information, quantum discord and entanglement (com- 
puted over the input two- mode state). 

As shown in the upper panels of Fig. [5J the QCB re- 
duction AQj is an increasing function of all the three 
correlation quantifiers for fixed input squeezing (/? = 0.1 
for the left panel and j3 — 0.9 for the right one). Note 
that, in each panel and for each quantifier, we plot three 
different curves corresponding to different values of the 
damping T = 0.9, 0.5 and 0.1. The monotonicity of 
the QCB reduction in all the correlation quantifiers sug- 
gests that the presence of correlations should definitely 
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FIG. 8: (Color online) Upper panels. QCB reduction AQj 
(with 7 = 0.999) as a function of the three correlation quan- 
tifiers X=I,D,E where I is the quantum mutual information 
(dotted red), D is the quantum discord (dashed blue) and 
E is the entanglement (solid black). The plots are for fixed 
squeezing: /3 — 0.1 for the left panel and ft = 0.9 for the 
right one. For each quantifier we plot three different curves 
corresponding to different values of the damping (from top 
to bottom F = 0.9,0.5 and 0.1). Each curve is generated by 
varying the input energy N between and 5 photons. Middle 
panels. Density plots of the QCB reduction AQ^ as a function 
of the input discord and entanglement. The plots are for fixed 
damping: T = 0.2 in the left panel and T = 0.8 in the right 
one. In each panel, the density plot is generated by varying 
the squeezing < /3 < 1 and the energy < N < 5. Lower 
panels. Entanglement (left) and discord (right) as a function 
of the quantum mutual information. Plots are generated by 
taking a random sample of 10 4 two-mode squeezed thermal 
states, i.e., random values of N and ft with 7 = 7. 



be considered as a resource for loss detection, whether 
these correlations are classical or genuinely quantum, i.e., 
those quantified by entanglement. In other words, em- 
ploying the input squeezing in the form of correlations 
is always beneficial for loss detection when we consider 
squeezed thermal states as input sources. The impor- 
tance of correlations is confirmed by the plots in the mid- 
dle panels. Here we consider again the QCB reduction 
AQ7 = Qi(N,/3) - Q 2 (-/V,/3,7) for 7 = 0.999. Then, by 
varying input squeezing j3 and energy TV, we study AQ^ 
as function of both discord and entanglement (damping 
is r = 0.2 in the left panel, and T = 0.8 in the right 
one). These plots show how the QCB reduction is ap- 
proximately an increasing function of both discord and 
entanglement. Finally, in the lower panels of Fig. [8j we 
also show how entanglement (left) and discord (right) 
are increasing functions of the quantum mutual informa- 
tion with good approximation (these plots are generated 
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by choosing a random sample of 10 4 two-mode squeezed 
thermal states). 

VIII. CONCLUSIONS 

In this paper we have considered the quantum discrimi- 
nation of lossy channels. In particular, we have focused to 
the case when one of the two channels is the identity, i.e., 
the problem of discriminating the presence of a damping 
process from its absence (loss detection). For this kind 
of discriminination we have considered thermal probes as 
input, i.e., single- and two-mode squeezed thermal states. 
The performance of the channel discrimination has been 
quantified using the QCB, computed over the two pos- 
sible states at the output of the unknown channel for 
a given input state. Finding the optimal input state p 
which minimizes this bound gives automatically the op- 
timal multi-copy state p<g> p® ■ ■ ■ when we consider many 
accesses to the unknown channel (under the assumption 
of single-copy constraints). In this scenario, we have 
fixed the total energy of the input state and optimized 
the discrimination (detection of loss) over the class of 
single- and two-mode squeezed thermal states. We have 
found numerically that the optimal states are pure, thus 
corresponding to single- and two-mode squeezed vacuum 
states. Furthermore, we have determined the conditions 
where the two-mode state outperforms the single-mode 
counterpart. This happens when the energy exceeds a 
certain threshold, which becomes zero for suitably low 
values of the transmissivity (i.e., high values of damp- 
ing). 

It is worth noticing that our approach (where we fix 
the total energy of probing and reference modes) also 
gives a sufficient condition for the problem where only 
the probing energy is fixed. In fact, if a two- mode state 
outperforms a single-mode state above a certain thresh- 
old value Nth of the total energy, this also happens when 
just the energy of the probing mode is above that value 
Nth- This is a trivial consequence of the fact that the 
total energy is bigger than the probing energy for two- 
mode states (A^ > iV2 robc ) while the two quantities are 



the same for sing le-mode states (JVi = iVf obe ). Thus, 

jyprobe = ^probe > ^ c£m be written ag JV 2 > AT X > 

Nth which is a stronger condition than N2 — Ni > N t h, 
since the QCB is decreasing in the total energy, as one 
can see from Eqs. and 

In our investigation we have then considered the prob- 
lem of loss detection in more realistic conditions, where 
it is unlikely to have pure squeezing. In this case, we 
have studied the optimal state for fixed total energy 
and squeezing, i.e., by fixing all the relevant resources 
needed to create the input state. Under these constraints, 
we have shown that a two-mode squeezed thermal state 
which conveys all the thermal photons in the dissipativc 
channel is the optimal thermal probe. In addition, this 
result is robust against fluctuations, i.e., it holds approx- 
imately also when the thermal photons are distributed 
in a more balanced way between the probe mode (sent 
through the dissipative channel) and the reference mode 
(bypassing the channel). 

Finally we have closely investigated the role of corre- 
lations in our problem of loss detection. We have found 
that, for fixed input squeezing, the reduction of the QCB 
is an increasing function of several correlation quantifiers, 
such as the quantum entanglement, the quantum discord 
and the quantum mutual information. We then verify 
that employing the input squeezing in the form of corre- 
lations (quantum or classical) is always beneficial for the 
detection of loss by means of thermal probes. 

The results of our paper provides new elements in the 
field of quantum channel discrimination and can be ap- 
plied to a wide range applications, including the char- 
acterization of absorbing materials. In particular, they 
are relevant in all the situations where the physical con- 
straints regard the creation of the input resources rather 
than the channel to be discriminated. 
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